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We report an ab initio study of the electronic properties of surface dangling-bond (SDB) states 
in hydrogen-terminated Si and Ge nanowires with diameters between 1 and 2 nm, Ge/Si nanowirc 
heterostructures, and Si and Ge (111) surfaces. We find that the charge transition levels e(+/— ) of 
SDB states behave as a common energy reference among Si and Ge wires and Si/Ge heterostructures, 
at 4.3 ± 0.1 eV below the vacuum level. Calculations of e(+/— ) for isolated atoms indicate that 
this nearly constant value is a periodic-table atomic property. 
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When a junction of two distinct, undoped semiconduc- 
tor materials (a heterojunction) is formed, one can define 
the band-lineup using the valence- and conduction-band 
edge energies on each side of the junction. The measure- 
ment or calculation of such lineup can be made either 
directly from the band discontinuities at the interface, 
or indirectly through a common bulk energy reference on 
both sides of the interface. From the theoretical side, sev- 
eral such common energy references have been proposed. 
Tersoff suggests [l| that an effective midgap energy Eb 
(which can be calculated from the bulk electronic struc- 
ture) can be used to predict lineups within an accuracy 
of 0.2 eV. Electronic states of transition metal impuri- 
ties 0-01 and charge transition levels £(+/—) of hydro- 
gen interstitial impurities 0, Q have also been shown to 
work as common energy references. 

In the present work, we report an ab initio study of 
the electronic properties of surface dangling bond (SDB) 
states in hydrogen-terminated Si and Ge thin nanowires 
(with diameters between 1 and 2 nm), Ge/Si nanowire 
heterostructures, and Si and Ge (111) surfaces. Here 
we define the SDB as the defect resulting from the in- 
complete passivation of a Si or Ge surface atom, due to 
a missing hydrogen. These defects are observed in hy- 
drogenated Si surfaces, and can be used as "gates" in 
molecular electronics devices 0]. They are also expected 
to occur in the recently synthesized ultra-thin hydrogen- 
terminated Si nanowires Q. We find that the charge 
transition levels e(+/— ) of SDB states are deep in the 
bandgap for the Si surface and wires, shallow for the Ge 
wires, and resonant in the valence band for the Ge (111) 
surface. We also find that the SDB e(+/— ) levels behave 
as a common energy reference level among Si and Ge 
wires and Ge/Si heterostructures, at 4.3 ± 0.1 eV below 
the vacuum level. Calculations of e(+/—) for Si and Ge 
isolated atoms, as well as for group III and V atoms, in- 
dicate that this constant value is a periodic-table atomic 
property. 

Our calculations are performed in the framework of 
Kohn-Sham density functional theory (DFT) [jj]. within 
the generalized-gradient approximation (GGA) [10( and 



norm-conserving pseudopotentials [111 ] in the Kleinman- 
Bylander factorized form [12J . We use the LC AO method 
implemented in the SIESTA code with a double- 
zeta basis set plus polarization orbitals. Total-energy 
differences (forces) are converged to within 10 meV/atom 
(40 meV/A) with respect to calculational parameters, 
that are the same as those of our previous work [Til ]. 

We consider nanowires based on the diamond bulk 
phase, oriented along the (110), (111), and (112) 
directions, with one SDB per unit cell. The H- 
terminated nanowire structures are obtained by choosing 
the nanowire axis along a given crystalline direction in 
the bulk lattice, and including all atoms that fall within 
a specific distance from the axis. The lowest-coordinated 
surface atoms are removed, and the dangling bonds of 
the remaining surface atoms are saturated with hydro- 
gen. Depending on the nanowire diameter and orienta- 
tion, three types of H passivation are possible, related 
to whether the surface Si (Ge) is bonded to one, two, or 
three H atoms, forming SiH (GeH), SiH 2 (GeH 2 ), or SiH 3 
(GeHs) radicals on the surface, respectively. 

Besides "pure" Si and Ge nanowires, we also consider: 
coaxial Ge/Si nanowires with a Ge core and a Si over- 
layer; a nanowire with a periodic alternation of a few 
Ge and Si layers along the wire axis, forming a Ge/Si 
nanowire heterostructure; and three-layer-thin Si and Ge 
films, hydrogen-terminated on both sides, representing 
H-terminated (111) surfaces. The structures of the coax- 
ial Ge/Si nanowires are shown in Fig. [TJ The structures 
of the pure Si and Ge nanowires have the same initial 
geometry as those in Fig. [TJ with all Ge and Si atoms 
replaced by either silicon or germanium. The geometry 
of the Ge/Si nanowire heterostructure is shown in Fig. [21 
with a periodic unit along the axis formed by the alter- 
nation of four layers of Ge atoms and four layers of Si 
atoms. The atomic positions and lattice parameter of 
all nanowires and (111) surfaces were relaxed within the 
tolerances indicate above. 

We now address the electronic structure of nanowires 
with one SDB defect per unit cell, obtained by removing a 
hydrogen atom from the nanowire surface and performing 
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FIG. 1: (Color online) Cross section of several Ge/Si coaxial 
nanowire heterostructures. In (a), (d), and (e), nanowircs 
with the axis along the (110) direction and with diameters 
(D) of 23 A, 15 A and 12 A, respectively; in (b), nanowires 
with the axis along the (112) direction and D — 19 A; in (c), 
nanowires with the axis along the (111) direction and D — 17 
A. In all cases, the nanowire axis is perpendicular to the plane 
of the paper. Si, Ge, and H atoms are represented by green, 
yellow, and blue circles, respectively. 
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FIG. 2: (Color online) Ge/Si nanowire heterostructure with 
axis oriented along the (110) direction and D = 1.2 nm. Side 
(a) and top (b) views of the unit cell are shown. 



full ab initio structural optimizations. We ensure that the 
unit cells we use are large enough to decouple neighboring 
defects. For the majority of the nanowires, we consider 
SDB's generated from GeH and SiH radicals, with the 
exception of wires with D — 1.4 (1.5) nm for which SiH2 
(GeH 2 ) radicals are considered [no SiH (GeH) radicals 
are present in this case]. Generally, we find that SDB 
Kohn-Sham eigenvalues are deep in the bandgap in the 
case of Si wires, and shallow (near the valence band edge) 
in the case of Ge wires. This is seen in Fig. [3l that shows 
the calculated band structures for the defective Si and 
Ge nanowires with D m 2 nm, oriented along the (110) 
direction. The unit cells are 7.9 A and 8.4 A long along 
the axis, for the Si and Ge wires, respectively. Both band 
structures in Fig. [3] are essentially identical to those of 



the perfect wires (not shown), except for the appearance 
of a half-occupied nearly-dispersionless band in the gap. 
The small dispersions (0.03 eV for Si and 0.09 eV for 
Ge) are characteristics of localized SDB states. These 
features of the SDB level are common to all ten Si and 
Ge "pure" nanowires we consider. 




FIG. 3: Electronic band structures of Ge (left side) and Si 
(right side) nanowires (D ~2 nm) with one surface dangling 
bond defect per unit cell. Both wires are oriented along the 
(110) direction. 

From Janak's theorem, Si = dE/drii [15J, and the 
nearly linear dependence of the SDB eigenvalue on its 
occupation number m, it follows that the eigenvalue of 
the neutral SDB state at the GGA, non-spin-polarized 
level, is a good approximation to the e(+/— ) charge 
transition level [16|. The latter is given by the value 
of the electronic chemical potential at which +1 and 
-1 charge states of a defect have equal formation en- 
ergies. In Fig. [5] we show a comparison between the 
SDB neutral eigenvalue and the £•(+/—) level, for the Si 
(D = 19 A) and Ge (D = 20 A) (112) wires. The straight 
lines for each charge state in Fig. [4] are Ej orm (SDB) = 
E q tot (SDB) - Etot (perfect) + m + qfi, where E q tot (SDB) 
is the supercell total energy for a nanowire with an SDB 
in charge state q, fx is the Fermi level, and Etot (perfect) 
is the total energy for the perfect wire. The chemical po- 
tential for hydrogen, /in, is obtained from a Si (Ge) bulk 
and SiH4 (GeH4) molecule calculations. The eigenvalue 
of the neutral SDB state is plotted as a dashed vertical 
line, and is within 0.05 eV from the e(+/— ) level for both 
Si and Ge. We obtained very similar results for Si and 
Ge wires oriented along the (110) direction. In the fol- 
lowing, we present strong numerical arguments showing 
that the e(+/— ) level of nanowire SDB's is a common 
reference level among Si, Ge, and Ge/Si nanowires. 

Let us first consider the electronic structure of the 
Ge/Si coaxial wires of Fig. Q] with an SDB state at the sur- 
face. Tabic |T] indicates the positions of the neutral SDB 
level relative to the valence band maximum, E& — E v , 
and relative to the vacuum level, Ed- The table shows 
that Ed — E v varies from a "shallow" value of 0.13 eV 
to a "deep" value of 0.50 eV, although the SDB's are, 
in all cases, at the surface Si atoms of the coaxial wires. 
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FIG. 4: Formation energies of SDBs in a (112) Si wire (di- 
ameter of 19 A) and a (112) Ge wire (diameter of 20 A). The 
vertical dashed line indicates the Kohn-Sham eigenvalue of 
the neutral SDB states. 



In contrast, Ed is a common energy reference for all the 
SiH-terminated wires (i.e., all in Table U except for the 
wire with D = 1.5 nm, which is SiH2-terminated). 



TABLE I: Axis orientation, diameter, calculated bandgap, 
and defect level position (Ed) [absolute and relative to the 
valence band maximum (E v )] of Ge/Si coaxial nanowire het- 
erostructures with one surface dangling bond per unit cell. 
Energies are given in eV. 



direction 


(110) 


(112) 


(111) 


(110) 


(110) 


Diameter 


23 A 


19 A 


17 A 


15 A 


12 A 


AEgap 


1.05 


1.40 


1.52 


1.35 


1.67 


Ed 


-4.23 


-4.25 


-4.23 


-4.43 


-4.24 


(Ed — E v ) 


0.13 


0.15 


0.29 


0.44 


0.50 



To further investigate the nature of the SDB states, 
we performed a calculation for a Ge (110) nanowire with 
D = 1.2 nm in which a GeH radical at the surface was 
replaced by a non-passivated Si atom, forming a Ge wire 
with an SDB at a substitutional Si atom at the surface. 
The band structures of this wire and of the corresponding 
one for a Ge wire with a "standard" Ge SDB are very 
similar to each other and to the one in Fig. (Ha) . The 
SDB state at the Ge atom is at E v + 0.27 eV, while the 
SDB state at the Si atom is at E v + 0.40 eV. The same 
is observed for a Si (110) nanowire with D = 1.2 nm in 
which a SiH radical at the surface is replaced by a non- 
passivated Ge atom: in this case, the SDB state at the Ge 
atom is at E v + 0.70 eV and the SDB state at the Si atom 



is at E v + 0.84 eV. These two cases demonstrate that the 
positions of the band edges relative to the SDB level are 
weakly dependent (within 0.15 eV) on the atomic species 
where the SDB is located, but strongly dependent on the 
electronic structure of the wire as a whole. 

Our third numerical argument in favor of the universal- 
ity of the SDB state is based on a Ge/Si heterostructure 
nanowire, shown in Fig. [2j We performed two separate 
calculations for SDB's in this wire. In the first one, a hy- 
drogen atom is removed from a Si surface atom. In the 
second one, we removed a hydrogen atom from a Ge sur- 
face atom. In both cases we placed the dangling bond as 
far as possible from the Ge/Si interface. Again, a strik- 
ing resemblance is observed in both band structures: the 
SDB level at the Ge atom is at E v + 0.51 eV while at 
the Si atom it is at E v + 0.59 eV. 
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FIG. 5: Calculated band edges often Si and Ge nanowires, rel- 
ative to the vacuum level. The top line indicates the conduc- 
tion band minima and the bottom line indicates the valence 
band maxima. The dashed line corresponds to the eigenvalue 
of the SDB state in each case. The orientation and diameter 
of each wire is also indicated. 

Given the above arguments, the e(+/—) level of the 
nanowire SDB's (or, likewise, its Kohn-Sham eigenvalue 
in GGA or LDA approximations) should be a common 
reference level among Si, Ge, and Ge/Si nanowires, and 
could be used, for instance, to predict band-lineups in 
nanowire junctions [l7l. flij and heterostructures. This is 
demonstrated in Fig. [3] where the calculated band edges 
(relative to the vacuum level [l9j]) of all pure Si and Ge 
wires are shown together with the eigenvalues of the SDB 
states. Figure [5] shows that the SDB e(+/— ) is a com- 
mon energy level among the nanowires. The similarity of 
the £■(+/—) is even stronger if we remove from the com- 
parison the SiH2- and GeH2-terminated nanowires (the 
ones with diameters of 1.4 and 1.5 nm, respectively). By 
doing that, all £■(+/—) values fall in a narrow range from 
-4.40 eV to -4.32 eV. 

We also performed calculations for SDB's on H- 
terminated Si and Ge (111) surfaces. We obtained 
e(+/— ) = —4.41 eV (relative to the vacuum level) for 
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the SDB on a Si surface, consistent with our results for 
the Si nanowires. For the SDB on a Ge surface, we find 
a resonant £(+/—) within the uppermost valence band, 
consistent with the fact that the top of the valence band 
E v = —4.26 eV is above the e(+/— ) nanowire values. In 
the case of the calculated Ge nanowires, Fig.[5]shows that 
E v is always below -4.4 eV (and therefore below e(+/— )), 
with a a clear trend of increasing E v with diameter D. 
This indicates that the E v values in Ge nanowires are 
downshifted in energy relative to the surface E v value 
due to quantum confinement. Therefore, for some di- 
ameter D (larger than the ones considered in this work) 
there will be a transition between localized and resonant 
e(+/— ), with large-diameter Ge wires presenting a reso- 
nant £(+/—}■ This is consistent with recent experimen- 
tal results [20] showing that large-diameter Si-covered Ge 
wires are intrinsically p-type. 

Why should the SDB's £(+/—) be a common en- 
ergy reference? Previously, it has been proposed that 
both in the case of transition- metal impurities and 
of interstitial- hydrogen impurities the common en- 
ergy reference is tied to an average dangling-bond energy 
which would be independent of the host material. This 
proposal has never been directly verified from first prin- 
ciples calculations, possibly because isolated dangling 
bonds are too reactive in bulk point defects. We argue 
that a surface dangling bond resulting from the removal 
of a hydrogen atom from a surface GcH or SiH radical 
should be a good realization of an ideal, unreconstructed 
dangling bond, as the surface Si or Ge atom with an SDB 
is bound to three saturated bulk-like atoms, with bond 
angles (9) and lengths (d) that are very similar to the 
defectless wire values (Ad < 0.03 A; A9 < 2°, 6°, and 
1° for Si, Ge and Ge/Si wires, respectively). Such small 
relaxations are indicative of small coupling between the 
SDB and the underlying bulk states. 

To gain an even deeper insight into the physical ori- 
gin of the SDB e(+/— ) universality, let us consider 
this quantity in the limit of isolated atoms. We com- 
puted the atomic-limit e(+/— ) using the same DFT ap- 
proximations as in the nanowire calculations. In this 
limit, e(+/— ) is defined as the average between the first- 
ionization potential and the electron affinity, computed 
for the atomic ground state. We obtained £•(+/—) = 
-4.39 eV for Si and -4.17 eV for Ge. This suggests 
that the universality is an atomic property. From simi- 
lar atomic calculations for group III and V neighbors to 
Si and Ge, namely, Al, Ga, P, and As, we computed the 
atomic-limit average e(+/—) = [e(+/— )v]/2 
of the III-V compounds A1P, GaP, AlAs, and GaAs. We 
obtained e = -4.37 eV (A1P), -4.34 eV (GaP), -4.22 eV 
(AlAs), and —4.20 eV (GaAs), indicating that the appli- 
cability of the SDB common energy reference should go 
beyond Si and Ge. 



In conclusion, ab initio calculations for hydrogen- 
terminated Si and Ge nanowires, nanowire heterostruc- 
tures, and Si and Ge (111) surfaces, indicate that the 
charge transition levels e(+/— ) of surface dangling bond 
states behave as a common energy reference among Si 
and Ge wires and Si/Ge heterostructures, at 4.3 ± 0.1 
eV below the vacuum level. Calculations of £(+/—) for 
isolated atoms indicate that this nearly constant value is 
a periodic-table atomic property. 
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